Correlation energy contribution to nuclear masses 
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The ground state correlation energies associated with collective surface and pairing vibrations 
are calculated for Pb- and Ca-isotopes. It is shown that this contribution, when added to those 
predicted by one of the most accurate modern nuclear mass formula (HFBCS MSk7 mass formula) , 
reduces the associated rms error by an important factor, making mean field theory, once its time 
dependence is taken into account, a quantitative predictive tool for nuclear masses. 

PACS numbers: 



Many of the nuclei in atoms we find in nature are pro- 
cessed in the stars. In particular, heavy elements far 
beyond 56 Fe are formed through a chain of neutron cap- 
ture reactions with subsequent beta decay. The source 
of the rather low energy neutrons are reactions like (a,n) 
on 22 Ni and 13 C where the alphas are highly abundant 
during He burning in the stars. The small neutron abun- 
dance then causes a particular path to be followed which 
lies mainly along the valley of beta stability as the cap- 
ture process follows a time scale controlled by the va- 
riety of /3-decay processes. The end point of this slow 
processes (s-process) is above 209 Bi via a-decay. There 
is clear evidence for sharp peaks in the abundance curve 
near neutron closed shell configurations which commonly 
result from the existence of slow neutron capture pro- 
cesses. However, the abundance curve displays a number 
of additional peaks at around 8-12 mass units prior to 
the s-process abundance peaks. This indicates that other 
capture processes where neutron capture is involved have 
to take place to produce the corresponding elements as 
well as elements beyond 2Q9 Bi. This is called the rapid 
or r-process. It is known that the r-process path evolves 
along the neutron drip line region where the bind- 
ing energy of the least bound neutron is quite small (few 
hundred KeV) in the presence of a large neutron flux 
Q ■ The need for a formula able to predict known masses 
with an accuracy of this order of magnitude seems quite 
natural. On the other there still exist many uncertain- 
ties in extrapolating our present-day knowledge of for- 
mulas which predict nuclear masses, towards the region 
of the drip line where weakly bound nuclei are found. 
In fact, the best account of the experimental data based 
on mean-field theory provides a fitting to the 2135 mea- 
sured masses with N,Z > 8 with a rms error of 0.674 MeV 
This has been achieved by means of Hartree-Fock- 
Bardeen-Cooper-Schrieffer (HFBCS) calculations which 
employ a Skyrme-type zero-range effective force in the 
mean field channel, supplemented by a zero-range pair- 
ing interaction. The 14 parameters set is named BSk2. 
As a reference point for our work, we have considered a 
parameter set of almost equal quality, denoted by MSk7, 



where the rms error is 0.738 MeV 4]. Although this ac- 
curacy is remarkable, one is still not satisfied. In fact, 
while it takes tens of MeV to remove two nucleons from 
a nucleus lying along the valley of stability, this quan- 
tity becomes as low as 300 keV in halo nuclei like 11 Li 
5] . The need for a mass formula at least a factor of two 
more accurate than the MSk7 mass formula is evident. 

Nuclei display both single-particle and collective de- 
grees of freedom. Consequently, the corresponding 
ground states and associated nuclear masses reflect the 
zero point fluctuations (ZPF) associated with these 
modes. While mean field theory includes fluctuations 
associated with quasiparticles, it is only time depen- 
dent mean field theory which takes into account the zero 
point fluctuations associated with collective modes. De- 
spite this, the information in the literature concerning 
ground state correlation energies associated with col- 
lective modes is rather scarce. A recent suggestion to 
consider their effect was put forward by G.F. Bertsch 
and K. Hagino 6J. Later, realistic calculations for the 
quadrupole degree of freedom were performed for light 
nuclei and for a few selected isotopes within the so- 
called Generator Coordinate Method (GCM) 0. 

In the present paper we work out, making use of 
Random Phase Approximation (RPA) calculations, the 
ground state correlation energies associated with both 
surface (quadrupole and octupole modes) and pair- 
ing vibrations for the Ca- and Pb-isotopes. Because 
pairing vibrations [jj have a collective character only 
around closed shell nuclei (being essentially pure two- 
quasiparticle states lying on top of twice the pairing gap 
in superfluid systems, cf. also |Ty|), one expects the asso- 
ciated ZPF to lead to important corrections to the mass 
formula 3 J . This in keeping with the fact that the largest 
deviation from experiment found in this mass formula are 
observed in closed shell systems. It will be concluded that 
ZPF are important, in particular those associated with 
pairing vibrations in closed shell nuclei, in reducing the 
rms of the Skyrmc- HFBCS mass formulas. 

To derive the particle-hole RPA equations use can 
be made of the quasi-boson approximation where the 
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particle- hole operators a^a^ and (a^ai) + are replaced 
by the boson operators and T ki . Because collective 
vibration can be viewed as correlated particle-hole exci- 
tations, the corresponding boson creation operator can 
be written as 



and 



r+(n) = £pq£(n)r+ + *S(n)r w ) 



(i) 



where X]^ and Y k " are the forwardsgoing and backwards- 
going amplitudes fullfilling the normalization condition 



[T a (n),Ti(n)} = J2m i \ 2 -\Y k a l \ 2 ) = l. 



(2) 



hi 



The equations which determine the frequencies of the 
vibrational modes of quantum number a (with progres- 
sively high energy n = 1,2,...) are obtained from the 
relation 



[H,T+(n)] = tuj a (n)T+(n). 



(3) 



In the above equation H is the total hamiltonian, sum 
of a single-particle and a two-body interaction term. In 
the present calculations we have used self-consistently the 
Skyrme plus pairing interaction given by the MSk7 pa- 
rameter set. The RPA ground state energy is given by 
(cf. e.g. [HI) 

E RPA =E HF -(2X+l)J2 E \ Y ^ n ^ > ( 4 ) 

a,n hi 

in keeping with the fact that the amplitudes Y ki (n) are di- 
rectly related to the ground state correlations induced by 
the corresponding vibrational modes. The second term 
of the r.h.s. is called correlation energy. 

It is well known that open shell nuclei display a finite 
BCS gap A, as a result of the pairing correlations act- 
ing between pairs of nucleons moving in time reversal 
states close to the Fermi energy. In closed shell nuclei, 
pairing correlations are not strong enough to overcome 
the large single-particle gap, and A = 0. Nevertheless, 
pairing plays an important role in determining the (low- 
lying) structure of these nuclei. In particular, it leads 
to pairing vibrations 0, which change in two the num- 
ber of nucleons, that is vibrations which can be viewed 
as two correlated particles (pair addition modes) or two 
correlated holes (pair removal modes). An example of 
these modes is provided by the ground state of 210 Pb and 
of 206 Pb which can be viewed as the pair addition and 
the pair subtraction modes of 208 Pb respectively. These 
modes are specifically probed through two particle trans- 
fer reactions [T^ . The pair addition and pair subtraction 
modes can be written as 



r+(n) = E^»r^ + E F i°W r - 



(5) 



r+(n)^x»r+ 



(6) 



where = [a£a£]o creates a pair of nucleons coupled 
to angular momentum zero in levels with energy larger 
than the Fermi energy (sk > £f), while Ti = [aiaj]o an- 
nihilates a pair of nucleons in occupied states (e^ < £p). 
Use is made of a hamiltonian H = H sp + H p , sum of a 
single-particle term and of a pairing force with constant 
matrix H p = —G^2jj,apat,ajaj, the commutation re- 
lation given in Eq.© leads to the dispersion relations 
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2e k - hw a (n) 
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G 2e k + hw r (n) 
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2ej + Huj a (n) ' 



2a — huj r (n) ' 



(7) 



(8) 



where ft — (2j + 1)/2 is the pair degeneracy of the single- 
particle orbital with total angular momentum j, while 



e 3 ~ £ 3 

X a k (n) 
and 
Xl{n) = 



■ ep. The amplitudes are 



A a (n)VT^ 
2e k - hu a (n) ' 



A, 



2ei — hw r {n) 



Y t a (n) 



Ykin) 



A a (n)yn7 

2ej + Hu) a (n) ' 



(9) 



2ek + huj r (n) ' 



(10) 



A a (n) and A r (n) being normalization constants deter- 
mined from the relation given in Eq. J2J). 

In Fig. n we show the dispersion relations given in 
Eqs.O and JSJ calculated for 208 Pb for both protons 
and neutrons (cf. also 9]), making use of the valence 
orbitals of this nucleus. The valence orbitals used in our 
calculations were determined with the help of a Woods- 
Saxon potential with standard parametrization [15J and 
the energies have been replaced with the experimental 
values whenever available. The results obtained using 
the standard Woods-Saxon levels coincide with those cal- 
culated making use of the experimental energies within 
2%. Making use of the fact that the sum of the pairing 
binding energies of 2m Pb and 210 Pb as well as in 206 Hg 
and 210 Po are rs 2 MeV (in this last case one has to take 
into account the Coulomb repulsion between the two pro- 
tons, cf. e.g. 0]), one obtains the values of 2.7 MeV and 
2.2 MeV for the neutron pair addition and pair removal 
energies, the corresponding values for the proton channel 
bein g 3. 5 MeV and 3.1 MeV respectively (cf. refs 0, 
and [lj])- The associated Y's amplitudes are displayed 
in Tables d and HD 

Inserting these results in Eq.QJ, one obtains the 
ground state correlation energy values -0.399 MeV (neu- 
trons) and -0.449 MeV (protons) respectively. In all cal- 
culations we have kept the contribution of only the low- 
est (n = 1) pair addition and pair subtraction modes, 
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FIG. 1: Monopole-pairing-vibrations dispersion relation for 
(a) neutrons and (b) protons for the nucleus 208 Pb 
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0.1837 


4si/ 2 
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0.0886 


lj'l5/2 ® ljl5/2 


-0.1530 


3^5/2 


0.0915 


0.1864 


2d 5/2 ® 239/2 
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3^1/2 


0.7853 


0.0839 


3P1/2 ® 2/ 5/2 


0.8431 


2/5/2 


0.4841 


0.1209 


3^1/2 ® 3p 3/2 


0.3756 


3^3/2 


0.2879 


0.0899 


2/5/2 ® 2/ 5/2 


-0.2614 


1*13/2 


0.3347 


0.1402 


2/5/2 ® 3p 3/2 
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2/7/2 


0.1856 


0.0914 


3^3/2 ® 3p 3/2 


-0.1186 


lft.9/2 


0.1461 


0.0839 


3^3/2 ® 2/7/2 


0.1422 



TABLE I: RPA wavefunctions of the neutron pair addition 
(a) and pair removal (r) modes of 20S Pb with multipolaritics 
and parity A"" = + , 2 + . 



in keeping with the fact that, as a rule, the n 7^ 1 
modes are much less collective. Pairing vibrations with 
multipolarity A 7^ 0, in particular quadrupole and hex- 
adecapole pairing vibrations, have also been identified 
around closed shell nuclei (cf. E2, 0, and refs 
therein) . In Table IIIII we display the contributions to 
the ground state energy (i.e. E^pa as defined in Eq.QJ) 
associated with the monopole, quadrupole and hexade- 
capole pair addition and pair removal modes for both 
neutrons and protons associated with 208 Pb, the summed 
contribution amounting to -1.981 MeV (« -1.196 MeV 
-0.785 MeV). 

In Table ITVl we collect the corresponding contribution 
for a number of Pb-isotopes. Because pairing vibra- 
tions are collective modes only around closed shell nuclei, 
where particles and holes can be clearly distinguished, 
becoming non-collective two-quasiparticle modes outside 
closed shells (cf. e.g. ref. Q) we have considered the 
contribution of neutron pairing vibrations only for the 
closed shell system (while the proton pairing vibration 
were taken into account for all isotopes). Also shown 
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TABLE II: RPA wavefunctions of the proton pair addition 
(a) and pair removal (r) modes of 208 Pb with multipolarities 
and parity A^ = + ,2+. 
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-0.399 


-0.449 


-0.609 


-0.244 


-0.189 


-0.092 



TABLE III: Ground state correlation energies, arizing from 
the neutron (n) and protons (p) monopole, quadrupole and 
hexadecapole pairing vibrations in 208 Pb. 





M4 Pb 




2,18 Pb 


2W Pb 


2i2 Pb 


p-h vibrations 


-2.793 


-2.709 


-2.237 


-2.801 


-3.173 


pairing vibrations 


-0.785 


-0.785 


-1.981 


-0.785 


-0.785 



TABLE IV: Ground state correlation energies for the Pb- 
isotopes. 





il> Ca 


42 Ca 


4i Ca 


4b Ca 


4H Ca 


p-h vibrations 


-0.886 


-1.418 


-1.606 


-1.391 


-0.547 


pairing vibrations 


-4.761 


-2.978 


-3.239 


-3.500 


-5.823 



TABLE V: Ground state correlation energies for the Ca- 
isotopes. 





a (MeV) 




ref. [4] 


ref. [4] + RPA 


Pb 


0.646 


0.543 


Ca 


1.200 


0.466 


a 


0.964 


0.505 



TABLE VI: Root mean square error associate with the HF- 
BCS MSk7 mass formula of ref. |4| and to this formula (with 
sligthly adjusted parameters) plus the correlation contribu- 
tions associated with surface and pairing vibrations calculated 

fa 2 +a 2 \ 1/2 

in the RPA. The quantity a = I Ca 2 pb I is shown in the 
last line. 

in table IIVI are the contribution to Erpa arising from 
the low-lying collective particle-hole vibrations calculated 
making use of the MSk7 interaction to determine the sin- 
gle particle states and the particle-hole correlated modes 
(cf. Eq.QJ). Quadrupole and octupole vibrations with 
energy < 7 MeV, and exausting > 2% of the non-energy 
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weighted sum rule were included in the calculation of 
Erpa- These conditions essentially select the lowest 
(one-two) states displaying correlated wavefunctions. 




40 42 44 46 48 
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FIG. 2: Difference between the calculated and the experimen- 
tal ground state energies. The curve (A) refers to the HFBCS 
MSk7 results of ref. j4|. The curve (B) was obtained adding 
p-h and pairing correlations to the curve (A) and renormal- 
izing the Skyrme parameters as indicated in the text. 

Similar calculations were repeated for the calcium iso- 
topes 40 ~ 48 Ca. In this case, the contribution of the pro- 
ton pairing vibrations calculated for the closed shell sys- 
tems 40 Ca and 48 Ca were linearly interpolated for the 
other isotopes, and no hexadecapole modes were con- 
sidered. In Table E| we show the corresponding results, 
together with the contribution of the particle-hole vibra- 
tional modes. When adding the results of Tables Hvl and 



El to the HFBCS MSk7 mass formula of Ref. [|, the 
parameters of the Skyrme interaction should be refitted 
in order to provide the best reproduction of experimen- 
tal masses. This should be done on a large sample of 
isotopes and it is beyond the purpose of the present pa- 
per. If we restrict ourselves to Ca- isotopes (Pb-isotopes), 
a slight renormalization by a factor 0.9964 (0.99945) of 
the Skyrme parameters ti,Wo is enough to shift the cal- 
culated value of the masses upwards by w 4.3 MeV (3.5 
MeV) and minimize the rms deviation, giving the results 
displayed in Table IVll and, for the Calcium isotopes, also 
shown in Fig. [21 Averaging the rms deviations associ- 
ated with Ca- and Pb-isotopes, leads to a value of 0.505 
MeV as compared to the value of 0.964 MeV obtained 
making use of the results of ref. 0. Although a global 
readjustment of the mean field parameters should be en- 
visaged, the fact that the locally extracted rms deviations 
have been reduced by a factor of w 2 can be considered 
meaningful. 

Similarly to what was observed in assessing the role 
played by zero point fluctuations associated with sur- 
face and pairing vibrations in connection with the nu- 
clear mean square radius 0, 0, 0] , the alignment of 
rotational nuclei and the pairing phase transition as 
a function of angular momentum |20l 1 2 1 L |22| . we con- 
clude that zero point fluctuations play an important role 
in providing the A-dependent contributions to nuclear 
masses needed to make mean field theory, including its 
time dependence, a quantitative predictive tool. 
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